Augmented collisional ionization via excited states in XUV cluster interactions 
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The impact of atomic excited states is investigated via a detailed model of laser-cluster interac- 
tions, which is applied to rare gas clusters in intense femtosecond pulses in the extreme ultraviolet 
(XUV). This demonstrates the potential for a two-step ionization process in laser-cluster interac- 
tions, with the resulting intermediate excited states allowing for the creation of high charge states 
and the rapid dissemination of laser pulse energy. The consequences of this excitation mechanism 
are demonstrated through simulations of recent experiments in argon clusters interacting with XUV 
radiation, in which this two-step process is shown to play a primary role; this is consistent with our 
hypothesis that XUV-cluster interactions provide a unique window into the role of excited atomic 
states due to the relative lack of photoionization and laser field-driven phenomena. Our analysis 
suggests that atomic excited states may play an important role in interactions of intense radiation 
with materials in a variety of wavelength regimes, including potential implications for proposed 
studies of single molecule imaging with intense X-rays. 



Over the last decade, coherent radiation sources have 
been developed that can probe light-matter interactions 
at ever smaller wavelengths and unprecedented intensi- 
ties [TJ |2] . Recent experiments have moved into the ex- 
treme ultraviolet (XUV), including several that have ex- 
plored intense XUV interactions with rare gas clusters 
[2 H] . The XUV regime near 30 nm is unique because 
the photon energy is too small for inner shell ionization 
of rare gas atoms, yet too large (at reported intensities) 
for any appreciable laser field-driven processes that dom- 
inate intense laser-cluster interactions at longer wave- 
lengths, such as collisional heating. We therefore propose 
that this regime presents a unique opportunity to isolate 
the influence of the internal electronic structure of atoms 
within clusters on intense radiation-cluster interactions. 

To date, models of cluster interactions with intense 
radiation in the near IR [5] or in the vacuum ultravio- 
let (VUV) [6HS] have not incorporated the effects of the 
atomic structure of the constituent atoms. As a result, 
the influence of the atomic structure in such interactions 
is unknown over a broad range of wavelengths, including 
in the XUV. 

In this Letter, we introduce a general model that ex- 
plicitly incorporates the effect of atomic and ionic excited 
states on collisional ionization, and apply it to a molec- 
ular dynamics code for rare gas cluster interactions with 
intense radiation. By including the process whereby a 
ground-state electron of any charge species can be inter- 
nally promoted to an excited state by a colliding electron, 
ionization is allowed to occur through a two-step process: 
excitation followed by ionization from the excited state. 
This requires two sequential collisions, each of less energy 
than is required by single step collisional ionization from 
the ground state (the conventional model used widely 
across wavelength regimes). This generally applicable 
two-step process, which we call "augmented collisional 
ionization" (ACI), is applied to the 32 nm experiment 
with Argon clusters [3]. ACI is found to dominate this 



interaction, and allows the system to access higher ionic 
charge states than are obtained via the single step colli- 
sional ionization process alone. 

The incorporation of ACI into our laser-cluster inter- 
action model provides the first explanation of the obser- 
vation of high charge states (Ar 4+) in the 32nm system. 
Further, we find that cluster atoms quickly become ex- 
cited and then ionized through ACI, demonstrating that 
the pulse energy is spread rapidly through the cluster. 
This work therefore provides a fundamental insight into 
the dynamics of XUV-cluster interactions, an area that 
has, to date, been the subject of only a small number 
of theoretical investigations [3 |4j [10]. Several of these 
focused on the electron emission spectra of the 32 nm 
- Argon interaction, though did not account for the ex- 
perimentally observed high charge state. In particular, a 
multistep ionization model was proposed that explained 
the electron spectra for all but the highest pulse inten- 
sities [3 , whereas Ziaja et al obtained good agreement 
with the emission spectra [10 using a Kinetic Boltzmann 
model including non-equilibrium and equilibrium dynam- 
ics. 

Through our model, we find that atomic excited states 
in collisional systems play a crucial role in understanding 
both the charge state spectrum and how quickly charge 
states are produced due to increased energy transfer be- 
tween electrons and ions. Though the XUV provides a 
unique regime in which to study these effects, our results 
are of more general application. As a result, this work 
opens the door to the study of the importance of ex- 
cited states in other material systems and in wavelength 
regimes ranging from the infrared to X-rays. Further, 
the more rapid energy transfer to ions through the ACI 
process could also play an important role in explosion 
dynamics, which may be an important consideration in 
proposed single molecule imaging studies with intense X- 
rays [TT] . 

In this work, the motion of the particles is calculated 
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via classical molecular dynamics. Quantum effects are in- 
cluded by determining transition cross-sections and using 
a Monte-Carlo scheme to determine when excitation or 
ionization occurs. In the case of ionization, a free elec- 
tron is added to the dynamics simulation. In the case 
of excitation, no new electron is added, but the parent 
atom/ion is set to be in an excited state, which deter- 
mines its future ionization potential and cross-section. 

The photoionization probability is calculated at each 
time- iteration for every atom/ion, and this is used to 
test for photoionization events. The photon flux is deter- 
mined via the intensity, which is modulated both by the 
time profile of the pulse and by photon absorption. The 
cross-section for the photoionization of neutral Argon was 
obtained from Ref. [12], while those of ionic Argon were 
obtained using Los Alamos Atomic Physics Codes [T3] . 

The atom-electron collision in the cluster environment 
is modeled as an isolated atom-electron collision in a 
background cluster potential. This allows the collisional 
cross-sections for both the ionization and excitation to be 
calculated using an isolated atom-electron model. Single- 
step collisional ionization cross-sections were calculated 
using the semi-empirical Lotz formula [14]. The Born 
plane-wave approximation [15] was used for both exci- 
tation and ionization from an excited state. Computed 
collisional excitation cross-sections were compared to ex- 
perimentally measurable excitations to metastable states; 
our results for neutral Argon closely match experimen- 
tal results for the metastable part of the 3p 6 — >• 3p 5 4s 
excitation ( 3 P and 3 P 2 ) [IS]. 

The excited states used in this work consist of a subset 
of all possible excited states; only single electron excited 
states were used. Of those, only the lowest eight states 
with I < 4 were implemented. This subset is the most 
important, as it contains the lowest energy states. In- 
cluding more single electron excited states adds to the 
total collisional cross-section, although states near the 
threshold require almost as much energy as ionization 
and are thus almost as infrequent. The energy of the 
excited states and the cross-sections were obtained using 
the Hartree-Fock implementation of the Cowan code [15] . 

We find that the inclusion of excited states in our 
model does not alter photonionization. This is because 
photoexcitation away from a resonance is rare as only 
states of | Al | = 1 are accessible. At the same time, the 
photon wavelength must also be approximately equal to 
the transition energy. In electron-atom collisions, in con- 
trast, the electron can transfer any amount of its orbital 
angular momentum and its available kinetic energy. 

The influence of the cluster is accounted for by deter- 
mining the potential V p at the atom due to all particles 
farther than the nearest neighbor distance, R (> 6.4 Bohr 
for Argon), and approximating it as a constant external 
potential in the vicinity of the atom. The potential at 
the atom is reinterpreted as the cluster-free threshold of 
the atom, similar to Ref. [8]. 



The constant potential approximation is justified by 
considering an atom at r = that is surrounded by two 
neighboring ions on the z-axis at ±R. The true back- 
ground potential at the atom is given by 
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where £ = cos(0), Pi is the Legendre polynomial of order 
/, Zi and Z2 are the charges of the neighboring ions. 

Both neighbors contribute a constant term, Z/R, to 
the Hamiltonian of the atom which merely shifts the 
threshold. Using first-order perturbation theory, the low- 
est order correction to the eigenvalue is, 

AE = -^±^(cP \q\M, (3) 

where q is the quadrupole operator. The state 0o is an 
eigenstate of the atomic Hamiltonian, and therefore the 
angular momentum operator, making the dipole term 
zero. The potential of the electrons within R of the atom 
is also removed from V p . The collisional system therefore 
closely approximates an isolated system. The screening 
effects of the electrons within R on the ionization po- 
tential are not easily modeled microscopically due to the 
highly collisional nature of the system. Since ACI would 
not be fundamentally changed by screening effects, they 
are not included here. 

During a collision, the energetically accessible states 
are determined by the kinetic energy of the impacting 
electron relative to the target atom's cluster-free thresh- 
old. This includes ionization if energetically permissi- 
ble. The total cross-section for all of the possible final 
states is used to determine if any excitation or ioniza- 
tion event will occur. If the electron's impact parameter, 
b(= \v x r|/M), i s within this total cross-section, then a 
state is chosen at random, weighted by its cross-section 
relative to the total cross-section. 

In our model, we have not included any mechanisms 
for the decay of excited states. While it is possible for 
the excited states to decay radiatively, the shortest lived 
states will be states of |A/| = 1. These decays are on 
the order of nanoseconds [T7J [18] making them unlikely 
to have a large effect on our results. 

Collisional de-excitation, which entails an impact elec- 
tron gaining energy by de-exciting the atom, is more 
probable though still unlikely to play a role. This is 
because the cross-section for ionization from an excited 
state is much larger than the de-excitation cross-section, 
favoring ionization provided the impact electron has the 
small amount of energy needed to ionize. Further, even 
if we were to allow collisional de-excitation to occur, it 
would be unlikely to significantly alter our results, since 
any impact electron that would gain energy would more 
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than likely use that extra energy to excite or ionize an- 
other atom, leading to the same results, on average. 

To determine the importance of ACI, we simulated the 
interaction of Argon clusters exposed to 32.8 nm radia- 
tion both with and without the ACI mechanism included. 
Conventional collisional ionization was included in all 
simulations. We used a pulse width of 25 fs, and varied 
the intensity from 5 x 10 13 to 1 x 10 14 W/cm 2 . Though 
the reported average cluster size in Ref. [3 was 80 atoms, 
a population of larger clusters is possible in the cluster 
jet. Thus we considered two sizes: Ar 80 , and the closed- 
shell icosahedral structure Ari 47 . The smaller closed- 
shell icosahedral, Ar 55 , may also be present but would 
likely have a much smaller contribution to the Ar 4+ sig- 
nal. For each set of parameters, we performed 10 4 sepa- 
rate simulations, and we ran each Argo simulation to 800 
fs and each Ari 47 simulation to 1 ps. 

In Table [TJ we list the average yield of Ar 4+ per cluster 
for each set of simulations. The number on the left was 
obtained from the inner-ionized Ar 4+ observed at the 
end of the simulations (i.e., not taking into account the 
electrons that are loosely bound to individual atoms). 
The number on the right was obtained by excluding any 
4+ ions that have at least one bound electron at the end 
of the simulation; this then gives a more accurate picture 
of what would be observed at the detector. 

Table [i] shows that Ar 4+ only emerges at the intensity 
quoted in Ref. [3] when the ACI mechanism is included. 
Though we see some inner ionized 4+ states without ACI 
at higher intensities, the yield is orders of magnitude less 
than when ACI is included, and there were no 4+ that 
were bare ions, meaning they would be detected as lower 
charge states. In contrast, when ACI was included, some 
4+ always remained without bound electrons and would 
thus be detected as 4+. While we see an increase in the 
number of 4+ in the larger versus smaller clusters, the 
influence of cluster size is still smaller than the increase 



Ar 80 


ACI 

Intensity~^-~^ 


No 


Yes 


5 x 10 13 W/cm" 


0/0 


0.0114/.0020 


7.5 x lO 1 ^ W/cm 2 


0.0002/0 


0.060/0.0045 


10 x 10 13 W/cm 2 


0.0008/0 


0.149/.0099 


Ari47 




No 


Yes 


5 x 10 r6 W/cm 2 
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0.0008/0 


0.151/0.0258 


10 x 10 13 W/cm 2 


0.0037/0 
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TABLE I. Average yield of Ar 4+ at different intensities for 
Argon cluster sizes of 80 and (closed-shell icosahedral) 147 
with (left) and without (right) two-step augmented collisional 
ionization (ACI) mechanism. The numbers on the right side 
were obtained by excluding any 4+ ions that have at least 
one bound electron at the end of the simulation. 



from the ACI ionization channel. Thus we find that, at 
all intensities, ACI is the single largest effect. 

To further illustrate the dominance of ACI, Figure [I] 
plots that fraction of each ionization channel that led 
to the creation of the various charge states. The left 
histogram is for Argo while the right is for Ari4 7 . Pho- 
toionization contributes only to the creation of Ar +1 and 
Ar +2 , as expected. Since photoionization plays a lesser 
relative role in the larger clusters, we see that clusters be- 
come more collisional as their size increases. Single step 
collisional ionization without excited states contributes 
significantly only to the creation of Ar +1 . ACI, however, 
dominates the creation of Ar +2 to Ar +4 , and is the ex- 
clusive mechanism for the creation of Ar 4+ . 




Ar— >Ar l+ Ar'^Ar 2 * Ar^Ar 3 * Ar^Ar 4 * Ar— >Ar 1+ Ar ljl ->Ar 2+ Ar^Ar 3 * Ar^Ar 4 * 



FIG. 1. (color online) The relative contributions of photoion- 
ization (diagonal lines), single step ionization (boxes) and ACI 
(solid) in the creation of each ion species of Argon by the in- 
teraction of (a) Arso clusters and (b) Ari47 clusters with 32.8 
nm radiation at an intensity of 5 x 10 13 W/cm 2 . 

Next, we examine the dynamics of the interaction to 
further explore the role of ACI and the prevalence of 
the excited states. Fig. [5] plots the relative (inner ion- 
ized) charge state population as a function of time during 
the interaction for Argo clusters both without (top) and 
with (bottom) the ACI mechanism included. Prior to the 
peak of the pulse (at 34.3 fs), atoms are the predominant 
species regardless of whether or not ACI is present. Soon 
after the pulse peak the Ar 1+ species is the most preva- 
lent. This lasts until about the end of the pulse at 68.6 
fs when the Ar 2+ becomes the most prevalent. 

With ACI, the higher charge states appear earlier and 
in greater abundance. The peak in the Ar 1+ population 
happens roughly 2 fs earlier with ACI, the Ar 2+ surpasses 
Ar 1+ around 10 fs earlier, and a sizable Ar 3+ population 
emerges 7 fs earlier and surpasses the atomic popula- 
tion around 12 fs earlier than without ACI. The final 
abundance of Ar 3+ is more than an order of magnitude 
larger when ACI is present. In both cases the Ar 2+ is the 
most abundant species at the end, but the ratios between 
species populations are very different. With ACI we also 
see the emergence of Ar 4+ , which appears at the tail of 
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the laser pulse but reaches a plateau around 100 fs later, 
at which time it surpasses the abundance of atoms. 
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FIG. 2. (color online) The charge species population relative 
to the total number of atoms and ions as a function of time for 
Arso clusters interacting with radiation of intensity 5xl0 13 
W/cm 2 , both without (a) and with (b) ACL The insert in 
(b) gives the relative excited charge species population as a 
function of time. 

The insert in Fig. [2] displays the fraction of each charge 
species that is in excited states as a function of time. The 
abundance of the excited states (after the pulse peak) 
closely follows the abundance of the internal charge states 
and a sizable fraction of each charge species is excited. 

ACI thus allows the cluster to reach high charges and 
to do so earlier than would otherwise be possible. This is 
a result of having a significant proportion of each charge 
species excited. Consequently, impact electrons do not 
require much kinetic energy to ionize the excited elec- 
trons as compared to single step ionization or the ex- 
citation step itself. For larger clusters, in which more 
electrons are trapped by the cluster potential, the peak 
abundance of excited states occurs even earlier than in 
small clusters at the same intensity (not shown). 

In summary, we have shown that the detection of the 
Ar 4+ in Ref. [3 can be explained by augmented colli- 
sional ionization (ACI) , a two-step process wherein colli- 
sional ionization occurs via intermediate atomic and ionic 



excited states. ACI was found to dominate the cluster- 
XUV interaction dynamics. In particular, it occurs much 
more frequently than direct collisional ionization from 
the ground state, which alone cannot reproduce the high 
charge state observed in experiment. While we have ap- 
plied our model to cluster-XUV interactions, ACI will 
likely play a role in any system that is highly collisional, 
including in a variety of finite condensed systems inter- 
acting with intense radiation in a wide range of wave- 
length regimes. 

We have also shown that ACI provides a route to high 
charge states in a shorter time, due to the efficacy of the 
energy transfer by ACI afforded by the lower energy gaps 
of excitation compared to direct ionization. The energy 
imparted to a small finite system will thus be highly sen- 
sitive to the ramp up of the pulse, which may have con- 
sequences for the direct imaging of single biomolecules 
with short, intense Xray pulses. Thus a good understand- 
ing of excited states of the photo-sensitive parts of the 
biomolecule (i.e., atoms and molecules that have largest 
photoionization cross-sections), and a sharp ramp up of 
the imaging pulse will be crucial. 
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